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Abstract 

The paper is devoted to the study of transmission probelm between two 
Herschel-Bulkley fluids with different viscosities, yield limits and power 
law index. 
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1 Introduction 

The rigid viscoplastic and incompressible fluid of Herschel-Bulkley has been 
scrutinized and studied by mathematicians, physicists and engineers as inten- 
sively as the Navier-Stokes. While this model describes adequately a large class 
of flows. It has been used to model the flow of metals, plastic solids and a 
variety of polymers like paints. Due to existence of yield limit, the model can 
capture phenomena connected with the development of discontinuous stresses. 
The literature concerning this topic is extensive; see e.g. [I], [6], [7], [8] and 
references therein. 

The purpose of this paper is to formulate and prove the existence of weak 
solutions for a class of boundary transmission problems for the Herschel-Bulkley 
fluid. To this aim we consider a transmission problem between two Herschel- 
Bulkley fluids with different viscosities, yield limits and power law index. We 
suppose that there is no-slip at the contact interface. Such problem can model 
the superposition of two paints in drawing process. 

The paper is organized as follows. In Section 2 we present the steady-state 
mechanical problem of transmission between two Herschel-Bulkley fluids. We 
introduce some notations and preliminaries. Moreover, we derive the variational 
formulation of the problem. In Section 3, we are interested in the existence of 
weak solutions. 
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2 Problem Statement 



We consider a mathematical problem modelling the steady-state transmission 
problem between two rigid, viscoplastic and incompressible Herschel-Bulkley 
fluids flow. The flrst fluid occupies a bounded domain ili C K" {n = 2. 3) with 
the boundary dfli of class C^. The second one occupies a bounded domain 
^2 C M" {n = 2,3) with the boundary 9il2 of class C^. We denote by f2 the 
domain fii U fl2 and we suppose that 



where ro,ri,r2 are measurable domains and meas (Ti) ,meas {T2) > O.The 
fluids are acted upon by given volume forces of densities fi , £2 respectively. 

We denote by S„ the space of symmetric tensors on M". We deflne the inner 
product and the Euclidean norm on M" and §„, respectively, by 

u • V = uivi V u, V e M" and a - t = aimTim V cr, r e S„. 



Here and below, the indices i and j run from 1 to n and Einstein's convention 
is used. We denote by ct the deviator of a = {(Jim) given by 



where 5 = {6im) denotes the identity tensor. 

Let 1 < p < 2. We consider the rate of deformation operator defined for 



We denote by n the unit outward normal vector on the boundary Tq oriented 
to the exterior of f2i and to the interior of Q2, see the figure below. For every 
vector field v € W^'''' (fii)" we also write v for its trace on d^li, i = 1,2. 

The steady-state transmission problem for the Herschel-Bulkley fluids is 
given by the following mechanical problem. 

Problem PI. Find the velocity fields Uj = {uu) : — > M" and the stress 
field C7j = {aum) '■ — > S„, i = 1, 2 such that 



dni = Fo U Fi and = Fq U F2, 



|u| = (u-u)5 VueM" and \a\ = {a ■ a)^ Vcr e 



^ = {d'lm) , ^lm = (Jlm 

n 



every u e VF^-p (f])" by 



D (u) = (Am (u)) , Dim (u) = - [Ul^m + U/.m) • 



U2 • VU2 = divcr2 + f2 in 02- 



ui • Vul = divcTi + fi in Q,\. 



(2.1) 

(2.2) 




(2.3) 




(2.4) 
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divui = in rJi. (2.5) 

divu2 = in rJa- (2.6) 

ui = on Ti. (2.7) 

U2 = on T2. (2.8) 

ui - U2 = on To. (2.9) 

cTi • n - (T2 • n = on To. (2.10) 



Here, the flows in the domains Qi, ^2 are given, respectively, by equations 
(|2.ip and (|2.2p where the densities is assumed equal to one for the two fluids. 
Equations (|2.3p and (|2.4p represent, respectively, the constitutive laws of the two 
Herschel-Bulkley fluids where fj,2 > and (72 > are the consistencies and 
yield limits of the two fluids, respectively, 1 < pi, P2 < 2 are the power law index 
of the two fluids, respectively. (|2.5p and (|2.6p represent the incompressibility 
conditions for the two fluids, respectively. (|2.7p and (|2.8p give the velocities on 
the boundaries Fi and r2, respectively. Finally, on the boundary part Fq, (|2.9p 
and (j2.10p represent the transmission condition for liquid-liquid interface. 

Existence of weak solutions for the flow of Herschel-Bulkley fluid was shown 
in 1969 for p > ^^^j for which the energy equality holds and higher differentia- 
bility techniques can be applied, in 1997 for p > and recently for p > 
using the Lipschitz truncation method. Moreover, some existence results has 
been obtained for the thermal flow in 2010 concerning the case p > see [5], 
HI) 0, [U, [5] and [3]. Up to now, there are only a few results concerning the 
regularity of weak solutions, especially in three-dimensional domains. 

We us consider the following function spaces 

V{n{) = {v e VF^'P' (rJi)" : divv in and V = on F,} , 1,2. 



y = {(vi, V2) e F {Vli) X V {VI2) : vi - V2 = on Fq} . 
V (rJi) , i = 1, 2 is a Banach space equipped with the norm 



and V becomes a Banach space for the following norm 



Vi,V2 U/ = VI ,,.0,-1 + V2 



v(ni) ll^2ll\/(02) • 

For the rest of this article, we will denote by c possibly different positive 
constants depending only on the data of the problem and denoting by p' the 
conjugate of p. 

Let us introduce the functional i = 1, 2 and the operator cj) defined by 
B^ (vi, V2, V3) = y" vi • Vv2 ■ vsdx, ^ ^ 1, 2. (2.11) 
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ct>:V^V', (ui, U2) ^ (t> (ui, U2) : V (vi, V2) £ V 
{4> (ui, U2) , (vi, V2))y,^y =f^iJ\D (ui)l^i-' Z? (ui) • (vi) dx+ 

fi^ J \D (U2)|^^"' D (U2) ■ D (V2) da:. (2.12) 

We begin by recalling the following, which gives some properties of the con- 
vective operators Bi. 



Lemma 2.1 1. Suppose that 

3n 



n + 2 



<p^<2, i = 1,2. (2.13) 



Then, 

Bi, i = 1,2 is trilinear, continuous on V {ili) x V {ili) x V {Vti) . Moreover, 
V(vi,V2,V3) e V(n^) X V{n,) X V{n,) we have 

(Vl, V2, V3) + (Vi, V3, V2) = 

(_l)*+iy(vi.n)(v2-V3)d7o, z = 1,2. f^.i^; 

To 

where d'y^ represents the superficial measure on the boundary part Tq. 

2. The operator (j) is hemi- continuous, strictly monotone, hounded and co- 
ercive on V. 

Proof. 1. The proof of the continuity of Sj, t = 1, 2 on y (f2j) x F (fi^) x V (fi^) 
is an immediate consequence of Holder's inequality and Sobolev embeddings, 
see [9]. 

Concerning the equality (|2.13l) it is enough to use an integration by parts 
and using the incompressibility condition (|2.5p . (j2.6|) . the boundary conditions 
(|2.7I) . (|2.8|) and the transmission condition (|2.9p . 

2. We can easily prove that the operator cf) can be written 

((/)(ui,U2) ,(vl,V2))^^,^^^ = {dJi (D (ui)) , £> (vi))^p'^^^^^„^„^^p^^j^^^„^„ + 

(dJ2 {D (U2)) , D (V2))^pi(f,^)„x„,ip,(,,^)„x„ , 

where the functional , i = 1 , 2 is defined by 

J, : (fl^)^" C §„ ^ R, cr ^J, (ff) = ^ f lal"' dx, i = l,2, 

Pi J 
Hi 

and d represents the Gateaux derivate. 
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Furthermore, it easy to check that the functional Ji, i = 1,2 is convex 
and Gateaux differentiable on L^' (fij)^^" . Thus, dJi is hemi-continuous and 
monotone. The Gateaux derivate of J, at any point a G LP' (ri^)"^" is given by 



This leads after some algebraic manipulations, that for i = 1,2 

> f^ij {\dJi{(Ti)\-\dJi{a2)\){\ai\-\a2\)dx. 

Then, if ai ^ (T2 we find 

(dJi(ai) -dJi(c72),CTi -a2)^^j(j^,j„x„^^p,(j^,j„x„ >0. 

Which means that the functional dJi, i = 1, 2 is strictly monotone. 
Consequently, operator ^ is hemi-continuous and strictly monotone on V. 
Moreover, we obtain from the definition of (p 

PI P2 

|(?!>(U1,U2),(V1,V2)>^,^^| < Ail MviQ,) Mvin^)+^J'2 l|U2||^(n,) \W2\\v(n^) 

(PI P2 \ 

Mil|ui||?(ni)+/^2l|u2||^n^)j ||(vi,V2)||y V(ui,U2), (ui,U2)gV. 

Then, 

Pi P2 

U{ui,U2)\\y, <Ml||ui||^(f^^j+M2l|U2||^(n^) V(U1,U2) G V. 

Hence, (f> is bounded on V. 

Now, we find using the generalized Korn inequality 

{(l){ui,U2),{ui,U2))y,^y > c(||ui||^'(f^^) + ||U2||^'(f^^)) . 

It follows that 



ll(ul>U2)||y l|ui|ly(ni) + ll"2||y(Q,) 

By passage to the limit when ||(ui,U2)||y — > +00, we find 

(^!>(ui,U2),(ui,U2))y,^y rf 1 cosPi 6* + cosP2 6» 

— — > c lim = +00. 

||(ui,U2)||y r^-+oo,ee[o,f] rcos6' + rcos^ 

This proves that the operator (f> is coercive. 
Which permits us de conclude the proof. ■ 
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Remark 2.2 In (QP, the right hand side has sense, since the injection 

W^-^t'^' (Fo)" ^ (Fo)" , * = 1,2 

is continuous. In particular the trace application 

7o : VFi'P' ^ (ro)" , 1 = 1.2 

is continuous. 

3n 

From now on, we take < K < 2, i = 1, 2. The use of Green's formula 

n + 2 

under the conditions (2.5)-(2.10) permits us to derive the following variational 
formulation of the mechanical problem (PI). 

Problem P^. For prescribed data (fi,f2) G V . Find (ui,U2) £ V satisfying 
the variational inequality 

Bi (ui,Ui,Vi - Ui) + B2 (U2,U2,V2 - U2) + 
((/>(U1,U2) , (Vi - Ui,V2 - \l2))v'y.V 

+gi j \D{wi)\dx~gi j \D{ui)\dx + g2 j \D {■V2)\dx - g2 j \D{u2)\dx 

r2i Oi ^2 

> J ii- iwi-ui)dx + J {2- iv2-U2)dx V(vi,V2)eF. (2.15) 

3 Main Result 

In this section we establish an existence result to the problems (Pv)- 
Theorem 3.1 The problem (Py) admits a solution (ui,U2) G V. 
The proof will be done in two steps. 

First step. Take an arbitrary element (wi,W2) G V and consider the 
auxiliary problem. 

Problem ^'(wi,w2)- Find (ui,U2) = (ui (wi,W2),U2 (wi,W2)) G V solution 
of the variational inequality 

Bi (wi,Ui,Vi - Ui) + B2 (W2,U2,V2 - U2) + 
((/>(U1,U2) , (Vi - Ui,V2 - U2))y,xy 

+.91 J \D{vi)\dx-gi J \D{ui)\dx + g2 J \D{v2)\dx-g2 J \D{u2)\dx 
f2i Oi ^^2 ^2 

>yfi-(vi-Ui)dx + y"f2-(v2-U2)da; V(vi,V2)gF. (3.1) 



Oi 

and it satisfies the estimate 



(ui,U2)||^, < i?. (3.2) 
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Lemma 3.2 The problem (P[-wi,w2)) ^'^■^ unique solution 

(Ui,U2) = (Ui (wi,W2) ,U2 (wi,W2)) £ V. 

Proof. Let us introduce the operator 

•^^Iwi^ws) (U1,U2) I > 0(wi,W2) (U1,U2) : V(vi,V2) G V 

(•^(wi.ws) (U1:U2) , (V1,V2)^^^^^ = Bi (wi,Ui,Vi) + B2 (W2,U2,V2) + 

(0(ui,U2),(vi,v2))v,xV'. (3.3) 
First, we get using lemma 2.1 

Bi (wi,Ui,Ui) +B2 (W2,U2,U2) = i y |ui|^(wi -11)^70- 



i y" |u2|^ (w2 • n) rf7o V(ui,U2)eF. 



To 

The fact that (wi, W2) , (ui, U2) 6 1^ imphes that Wi — W2 = and Ui — U2 = 
on To. Thus, 

(wi,W2) (U1'U2),(U1,U2)^^^^^ = (0(ui,U2),(ui,U2))^,^y 

V(ui,U2) e V. (3.4) 
Furthermore, we find for every (ui,U2) , (vi,V2) S V. 

(0(wi.W2) (V1,V2) - 0(wi,W2) (U1,U2) ,(vi -Ui,V2 -U2)^^^^^ 
= Bi (wi,Ui, Vi - Ui) + B2 (W2,U2, V2 - U2) - Bi (Wi, Vi, Vi - Ui) - 
B2 (W2, V2, V2 - U2) + (0 (Vi, V2) - (j> (Ui,U2) , (vi, Vg) - (Ui , U2)) y, ^ y 
V(vi,V2) , (Ui,U2) e V. 

This gives, keeping in mind lemma 2.1 



'/'(W1.W2) (V1,V2) (Ui,U2),(vi -Ui,V2 -U2)y^, 

= ^ y"|ui-vi|^(w2-n)d7o-i J \U2 - V2\^ {w I ■ n) d-f 



V'xV 




To To 
+ (0 (V1,V2) - (/>(U1,U2) , (Vi - Ui,V2 - U2))y,^y 
= ('?!>(vi,V2) - (/i(ui,U2) , (Vi - Ui,V2 - U2))y,xv 

V(vi,V2), (ui,U2) e K (3.5) 

Consequently, lemma 2.1 leads making use ([1]) and (13. 5|) that the operator 
is (/>(wi,w2) hemi-continuous, strictly monotone, bounded and coercive on V 
for every (wi,W2) G V. 



7 



Consider now the following functional 

j(vi,V2)=.gi / \D{vi)\dx + g2 I \D{v2)\dx (3.6) 



We can easily verify that the functional j is proper, convex and lower semi- 
continuous on V. 

The inequality p.ip can be rewritten using the operator 0^^^ and the 
functional j as follows 

(W1.W2) (U1>U2) , (Vl - Ui,V2 - U2)^^^^^ + J (vi,V2) ~j (ui,U2) 

> y"fi • (vi -ui)dx + y"f2 • (v2 -U2)da; V(vi,V2)eF. (3.7) 

Oi n2 

Consequently, the existence and uniqueness results from classical theories 
for inequalities with monotone operators and convex functionals, see [I]. 

Furthermore the estimate p.2p can be easily deduced by setting (vi,V2) = 
(0,0) as test function in inequality (13.11) . using lemma 2.1, Korn's inequality 
and some algebraic manipulations. ■ 

Second step. In order to obtain the solution of problem (-Pv) from that 
of problem /'(wi,w2)i ^^e the Schauder fixed point theorem, see [S]. To this 
aim we introduce the ball 

= {(wi,W2) eV: ||(wi,W2)||v. < i?}, (3.8) 

where R is the constant given by the estimate (13.21) . The ball K is convex 
and from the Rellich compactness theorem the ball is compact in L"^ (^1)" x 
(1^2)" ■ Let us built the mapping C : K — > K, as follows 

(wi,Vif2) I > C{-Wi,V/2) = (U1,U2) . 

To conclude the proof it is enough to verify the continuity of the mapping C 
when the ball K is provided by the topology of space (f^i)" x (^2)" ■ 
To do this, we consider (wi,W2), (wi,W2) G K and denoting by (ui,U2), 
(u'j, U2) G K the elements (ui, U2) = £ (wi, W2) and (u'^, Uj) = £ (w^, Wj) . 

Remembering that (ui,U2) and {u[,U2) are the solution of the problems 
below 

Bi (wi,Ui,Vi - Ui) + B2 (W2,U2,V2 - U2) + 
((/>(ui,U2),(vi -Ui,V2 - U2))y,^y 

+gi J \D{vi)\dx-gi J \D{ui)\dx + g2 J \D{v2)\dx~g2 J \D{u2)\dx 

fij^ O]^ 

> J {i-{vi~ui)dx + J {2-iv2-U2)dx W{vi,V2)eV, (3.9) 

Oi SI2 
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and 



Bi (w;,u;,vi - U[) + B2 (w(,,U^,V2 - u^) + 
('/'(u'l,U2),(vi -U;,V2 -U2))y,xy 



+51 J \D{vi)\dx~gi J \D{u[)\dx + g2 J \D{v2)\dx-g2 J \D{u'2)\dx 
r^i Oi O2 ^2 

> y fi • (vi - u'l) dx + y f2 ■ (V2 - u'2) dx V (vi, V2) e V, (3.10) 

Oi £72 

Now, choosing Vi = u[ and V2 = U2 as test function in inequality p.9p 
and Vi = Ui and V2 = U2 as test function in inequality p.lOp . It follows by 
subtracting the two obtained inequalities and using lemma 3.1, the transmission 
conditions, the definition of the space V and some calculations 

Bi (w'l - Wi, Ui, u'l - Ui) + B2 (W2 - W2, U2, - U2) + 

2 



^ / '"'^ ~ "^'^ ^"^'^ • n) ^7o - ^ J |U2 - U2|^ {^'2 ■ n) d7o 

To To 

+Mi y (u'i)r^~' D K) - |D (uor^-^ (ui)) . D K - ui) 

+M2 / (U'2)r'"' K) - |I? (U2)r^-' (U2)) • D K ~ U2) dx < 0. 

(3.11) 

Observe that for every x,y G M", 

(\xr^x-\yr'y)-{x-y)>c-^p^, l<p<2. (3.12) 
Then, inequality p.lip becomes 

'J i\D{u,)\ + \D{u[)\f-''^ ^'J {\D{u2)\ + \Diu',)\f-''^ 

ill ''2 

< c\Bi (w'l - Wi,Ui,u'i - Ui)| +c\B2 (W2 - W2,U2,U2 - U2)| (3.13) 

On the other hand, the application of Korn's and Holder's inequalities leads 
for i = 1, 2 to 

II - Uilly(o.) 



i\Diu,)\ + \D{u'^\y 



(3.14) 
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This yields, taking into account (13.21) . p.l3p and Holder's inequality 



IKu'i - Ui,U2 - U2)||y < 

Thus, Sobolev' embedding leads via the estimate p.2p to 

c||(wl - wi,wi - W2)||^^^^^^„^^^^^^^„ . (3.15) 

Hence, by virtue of Schauder's fixed point theorem, the mapping C admits 
a fixed point (ui, U2) ~ C (ui, U2) , which solves the problem (Pv)- 
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